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ABSTRACT 


This  report  comprises  three  related  papers  on  inferential 
procedures  for  the  V’eibull  or  extreme-value  distribution  based  on 
censored  samples.  In  Chapter  I  a  simple,  unbiased  estimator, 
based  on  a  censored  sample,  is  proposed  for  the  scale  parameter 
of  the  extreme-value  distribution.  The  exact  distribution  of  the 
estimator  is  determined  for  the  cases  in  which  only  the  first  two 
or  only  the  first  three  observations  are  available.  The  asymptotic 
distribution  is  derived,  and  an  approximate  distribution  for  small 
sample  size  is  also  provided.  Interval  estimation  for  the  scale 
parameter  is  developed  and  a  conservative  interval  estimate  for 
reliability  is  also  obtained. 

Chapter  II  represents  a  continuation  of  the  material  in 
Chapter  I  with  emphasis  on  combining  independent  lots.  A  study 
of  the  saving  in  experiment  time  with  censored  sampling  and  a 
numerical  example  are  also  presented. 

In  Chapter  III  tables  are  provided  for  obtaining  confidence 
limits  on  tie  parameters  or  reliability  based  on  the  maximum  like¬ 
lihood  estimators  for  selected  censoring  and  sample  sizes. 
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CHAPTER  I 


INFERENCES  BASED  ON  CENSORED  SAMPLING 
FROM  THE  WE I BULL  OR  EXTREME-VALUE  DISTRIBUTION 


1.  INTRODUCTION 

The  Weibull  distribution  denoted  by 

Fx(x)  =  1  -  exp  [~(x/a)S],  0  <  X  <  °°;  cs,  8  >  0, 

is  considered  in  this  paper.  The  variable  Y  =  In  X  follows  the 
extreme-value  distribution 

F  (y)  =  1  -  exp  [-exp  [y-u)/b]  ]  ,  -“  <  y  <  «>, 

B 

where  b  =  1/8  and  u  =  In  a,  The  variables  (X/a)  and  (Y-u)/b 
are  the  corresponding  reduced  variates  whose  distributions  are 
given  by  letting  a  =  6  =  1  and  b  =  1,  u  *  0,  respectively. 
Equivalent  procedures  can  be  developed  under  either  model,  but 
the  extreme-value  distribution  has  the  advantage  that  its  para¬ 
meters  appear  as  location  and  scale  parameters. 

Point  and  interval  estimation  procedures  are  in  general 
quite  complicated  for  these  models,  especially  under  censored 
sampling.  Maximum  likelihood  estimators  of  the  parameters  can  be 
determined  with  the  aid  of  a  computer  [4,5];  and,  inferential 
procedures  based  on  the  maximum  likelihood  estimators  have  been 
rather  extensively  developed  for  complete  samples  [16,17].  Theo¬ 
retically  these  methods  can  be  extended  to  censored  samples; 
however,  the  amount  of  computer  time  needed  in  order  to  determine 
the  appropriate  distributions  becomes  excessive.  Thus  other 
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estimation  and  hypothesis  testing  techniques  need  to  be  developed 
for  the  censored  sampling  case. 

For  point  estimation  one  common  approach  has  been  to  apply 
the  generalized  least  squares  method  or  some  related  method  to 
obtain  linear  estimators  of  the  location  and  scale  parameters  v 
and  b ;  see  for  example  [9,10,11].  This  approach  for  the  most  part 
requires  knowledge  of  the  variances  and  covariances  of  the  ordered 
observations  of  the  extreme-value  distribution,  which  are  avail¬ 
able  up  to  n  =  25  [13].  It  does  not  appear  convenient  to  deter¬ 
mine  tests  or  confidence  intervals  based  on  these  point  estimators 
for  sample  sizes  larger  than  25.  Some  simple  alternate  point  and 
interval  estimation  procedures  have  been  presented  in  [11,14]  for 
censored  sampling.  Also,  Johns  and  Lieberman  [8]  have  a  notable 
paper  concerned  with  determining  lower  bounds  for  reliability 
based  on  censored  samples.  An  attempt  is  made  in  this  paper  to 
develop  procedures  which  are  simple,  reasonably  good  and  widely 
applicable  without  the  necessity  of  generating  an  undue  number  of 
tables . 


2 .  INFERENCES  CONCERNING  b 
2.1.  UNBIASED  ESTIMATOR  OF  b 

Suppose  x. , . . . , x  denote  the  i  smallest  ordered  observations 
in  a  sample  of  size  n  from  the  Weibull  distribution.  Also 
yi,...yr,  where  y^  *  In  x^f  will  represent  the  r  smallest  observa¬ 
tions  in  a  sample  of  size  n  from  the  extreme-value  distribution. 
The  corresponding  reduced  observations  will  be  denoted  by 

zi  “  i and  wi  * 
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On  examining  a  table  of  coefficients  for  determining  best 
linear  unbiased  estimators  (BLUE'S) of  b,  one  sees  that  the  statistic 


r-1 


b  = 


■  '.I  «VV/nkr,n  -  T/kr,n 
1  =  1  ' 


is  an  appropriate  unbiased  estimator  of  b,  where 
r-1 

k  =  -(l/n)E  £  (W.-W  ).  The  statistic  b  is  the  BLUE  for  r=2, 

r,n  i=i  1  r 

and  it  is  somewhat  similar  to  the  BLUE  for  larger  r.  The  exact 

moments  of  the  W.  are  given  in  [18]  up  to  n  =  100,  and  k  can 
x  r  f  n 

be  calculated  easily  from  these  values  for  any  prescribed  com¬ 
bination  of  r  and  n.  For  illustration  purposes,  values  of  k^  ^ 
are  presented  in  Table  I  for  n  =  5,10,15,20,30,60,100  and  integer 
values  of  r=np>2  for  p  =  .  1(. 1)1.0.  Asymptotic  results  can  be 
utilized  for  larger  values  of  n.  If  r/n  p  as  n  °»,  expressions 
for  the  asymptotic  values,  say  k  ,  of  the  constants  have  been 

Jr 

derived,  and  numerical  values  of  k  are  also  presented  in  Table  I. 

P 

Values  of  the  asymptotic  efficiency  of  b  are  also  given  in  Table  I, 
and  these  indicate  that  b  is  a  simple,  relatively  efficient  point 
estimator  for  b  under  censored  sampling.  A  comparison  of  the 
■Variance  of  b  with  the  variance  of  the  best  linear  unbiased  estimator 
of  b  is  also  given  in  Table  II  for  some  small  sample  sizes  and 
various  censoring  fractions.  It  may  be  worth  noting  that  if 
mean  squared  error  rather  than  variance  is  used  as  a  goodness 
criterion,  then  a  value  of  c  can  be  found  such  that  the  MSE(cb)  is 
minimized.  This  value  of  c  is  given  by 

c  =  [1+var (b/b) ]_1  =  nk  /(I+nk  ). 

r ,  n  r ,  n 

A  comparison  of  the  mean  squared  error  of  cb  with  the  mean  squared 
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error  of  the  best  linear  invariant  estimator  (BLIE)  [9]  of  b  is  also 
given  in  Table  II. 


2.2.  SUMMARY  OF  DISTRIBUTION  RESULTS 

In  considering  similar  tests  concerning  b  relative  to  the 
nuisance  parameter  u,  one  observes  that  attention  can  be  re¬ 
stricted  to  functions  of  the  r  -  1  statistics  Y.  -  Y  =  ln(X./X  ); 

i  r  i  t* 

i  =  l,...,r  -  1.  The  statistic  b  is  suggested  for  use  since  it 
is  a  natural  type  of  function  of  these  statistics  to  consider,  and 
it  is  also  a  good  point  estimator  of  b  under  censored  sampling. 
Another  possible  indication  of  its  desirability  is  that,  although 
it  is  not  a  sufficient  statistic,  it  does  appear  as  a  quantity  in 
the  joint  density  of  the  Y^  -  Y  .  Also,  since  b/b  is  distributed 
independently  of  all  parameters,  tests  based  on  b  are  convenient 
to  express,  if  the  appropriate  percentage  points  are  available. 

A 

The  derivation  of  the  distribution  of  S  =  exp[-nk  b/b] 

r-1  r,n 

u  6 

=  n  (x./x  )  is  considered  in  section  5.1.  For  r  =  2,  the 
i=l  1  r 

distribution  of  S  is  shown  to  be 


F2(s)  =  ns/(s+n-l),  0  <  s  <  1. 


For  r  =  3, 


F3  (s) 


n!  f  1  1  (n  -  2)  (s  -  ] 

i  -  3 )  1  |  2n  "  2  (s  +  n  -  1)  '  2q  (s  +  n  -  1] 


s 

i371 


-  /q)  (2s  +  n  -  2  +  /g)!} 
(n  +  /qj  (2s  +  n  -  2  -  /q)  ^ 


where  q  =  (n  -  2)  -  4s. 
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The  exact  distribution  becomes  intractable  for  larger  r; 

however,  the  chi-square  distribution  can  be  used  to  provide  an 

extremely  good  approximation.  It  is  shown  is  section  5.2  that 

-2  In  S  =  2nk  b/b  is  distributed  approximately  as  a  chi-square 
r ,  n 

variable  with  2nk  degrees  of  freedom,  for  r/n  about  .5  or 

r  ,n  3 

less . 

The  asymptotic  distribution  of  T/b  =  -In  S  is  also  derived 
in  section  4.  It  is  shown  that  the  distribution  of 
Sn((T/b)  -  ^p^/°p  approaches  a  standard  normal  distribution  as 
n  -*■  »  and  r/n  -*■  p,  where 

oo 

un  =  I  (-'V1/  (i)  (i)!' 

p  i=l  p 

°P  =  P3/(l  -  P)Ap  -  Mp  +  2 Upp/Ap  +  2£).p  (-1)  1+1/i2i I  , 

2 

and  A  =  -In  (1  -  p) .  Numerical  values  of  k  =  -y  and  a  are 
P  P  P  P 

tabulated  in  Table  I  for  p  =  . 1 ( . 1 X  9 . 

2.3  TEST  OF  HYPOTHESIS  CONCERNING  b  AND  POWER  OF  THE  TEST 

Consider,  for  example,  the  test  of  Hq:  b  <  bQ  against  the 
alternative  HA:  b  >  bQ  at  the  a  significance  level.  Using  the 
chi-square  approximation  stated  in  section  2.2,  one  rejects  the 
hypothesis  Hq  if 

-2  l  (Y.  -  Y  )/b  >  xf  (2nk  ),  where  Prix2(v)  >  x.2(v)]  =  a. 

^  X  X  v  G*  A  f  X 1  '■« 

Linear  interpolation  can  be  used  for  non-integer  degrees  of 
freedom. 

%• 
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The  power  of  the  test  for  an  alternative  b  is 

Pr  [reject  hq]  =  Pt-2^  -  Yr)/bQ  >  X^nk^)] 

-  pl*2(2nkr,n>  >  <Vb)*«(2nkr..n)1 

Note  that  a  test  on  b  is  analogous  to  a  test  on  the  variance  of  a 

2  2 

normal  distribution,  since  in  that  case  (n  -  l)s  /a  is  distri¬ 
buted  as  a  chi-square  variable  with  n  -  1  degrees  of  freedom. 

Thus  material  developed  for  the  normal  case  can  be  applied  to 

this  case.  For  example,  the  sample  size  table  and  o.c.  curves  on 

2 

pages  299-303  of  [2]  are  applicable  by  simply  replacing  a  by  b 

and  n  by  2nk  +  1. 

2  r  ,n 

3.  INFERENCES  ON  THE  RELIABILITY 

The  problem  of  determining  a  test  for  £  =  a1/*5,  or  equivalently 
the  reliability,  R  =  exp[-(t/a)  ],  will  now  be  considered.  It  is 

[r_l 

l  +  (n  -  r  +  1)X  is  a 

i=l  r  -1 

complete,  sufficient  statistic  for  a  if  b  is  known,  and  that 
2r£/£  is  distributed  as  a  chi-square  variable  with  2r  degrees  of 
freedom.  Since  the  distribution  of  b  is  independent  of  a,  it 
follows  from  Basu's  Theorem  [1]  that  £  and  b  are  stochastically 
independent.  Thus  a  joint  confidence  region  can  be  determined 
for  b  and  R,  and  a  conservative  limit  for  R  can  be  obtained.  A 
similar  approach  has  been  followed  by  Mann  [9]  to  obtain  a  test 
for  R  based  on  Xr/X^  and  2r(./i. 
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In  terms  of  R, 


2r 5/c  =  2  (-In  R) [Z (X±/t) 1/b  +  (n-r+l)  (xr/t) 1/b] , 


and 


[ki  <  x2 

L  5  X“! 


(2r),  <  2"kr  nb/b  <  x^(2nkr  ^ 


=  (1  -  a1)  (1  -  aj  ”  “3)  * 


This  gives  the  joint  confidence  region  {R  >  R(b) ,  b  <  b  <  b}, 
where 

R(b)  =  exp(-xf  (2r)/2[J  (Xi/t.)1/b  +  (n-r+l)  (Xr/t) 1/b] ) , 


k-  2nkr,ni/xa3  (2nkr,n>’  and  b  *  Mr,nl/<L.2(Mr,n>' 


A  conservative  (1  -  a-^)  (1  -  ~  a^)  confidence  limit  for  R  is 

then  given  by  R  =  min  __  R(b)  .  It  is  shown  by  Mann  [9]  that 
”  b<b<b 

R(b)  is  a  monotonically  decreasing  function  of  b,  if  the  time  t 

r  i/r 

is  sufiiciently  small,  and  at  least  if  t  <  (  n  x. )  ,  in  which 

i=l 

case  R  =  R (b) .  This  would  ordinarily  be  the  situation  if  p  is 

not  too  small,  since  the  expected  coverage  between  any  two  order 

statistics  is  l/(n‘l).  It  is  also  clear  that  if  t  >  x ,  then  R(b) 

is  an  increasing  function  of  b,  and  R  =  R(b).  If  t  <  x  but  near 

xr,  then  R(b)  may  not  be  monotonic;  but  it  would  have  a  single 

minimum,  and  R(b)  would  approach  1  as  b  +  0,  and  R(b)  would  approach 
2 

exp{-x  (2r)/2n)  as  b  -*•  «.  Thus,  if  t  <  x_  and  R(b)  is  de- 

*  ft  i 

®1 

creasing  at  b,  then  R  =  R(b) .  If  P(b)  is  increasing  at  b,  then  a 
search  between  b  and  b  would  be  needed  to  determine  the  minimum. 
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4.  ASYMPTOTIC  DISTRIBUTION  OF  b 


Results  given  in  [3]  will  be  applied  to  determine  the  asymp- 

r- 1 

totic  distribution  of  -T  =  l  (Y.-Y  )/n,  where  r/n  ■*  p  as  n  ■+■  °°.  If 

i=l  1  r 

r 

we  follow  the  notation  of  [3],  -nT/b  =  F  In  X.  -  r  In  X 

„  1 
=  l  c.  h(X.)/  where  the  X.  are  ordered  exponential  variables. 
i=l  n  1  1 

Also,  oin  -  1.  i  -  1 . r-1,  crn  =  -r,  cin  =  o,  i  -  r  +  1 . n; 

F  (x)  =  1  -  exp(-x),  h(x)  =  In  x  =  H(x),  H(u)  -  h[F~1(u)]  =  In  [-In  (1-u)  ]  . 
Also 

i 

E(X  )  =  l  1/ (n-j+1) , 

j=l 

n  ~ 

a.  =  [1/ (n-j+1)]  l  c.  H*  (E  (X ■ ) ) 
j  n  .  .  in  l 

1=3 
r— 1 

=  (n-j+l)'1!  J  [E(Xi)]"1  -  r[E(Xr)]-1]. 


Then 


2 

a 


n 


(1/n)  l  c.  H[E(X.)3 
j  =  l  -|n  3 

r 

(1/n)  l  In  E(X. )  -  (r/n)  In  E(X  ), 
j  =  l  3  r 

U/n)  l  ^ 

j=l  3n 


and  the  distribution  of  /n((-T/b)  -  un)/a  approaches  a  standard 

2  2 

normal  distribution.  Furthermore,  u  m  and  a  -*•  o  ,  where 

n  p  n  p 

r/n  ■*  p  as  n  -*•  00 .  Let 


J (x)  =1,  x  <  r/  (n+1) 
=  0,  otherwise, 
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and  let  a,  =  p,  A  =  F  1 (p)  =  -ln(l-p).  Then 

***  ir 


J  (u)  H  (u)  du  -  a,  h  ( A  ) 

Ip 


In [-In  (1-u) ] du  -  p  In  A 


=  l  (-X  )V(i)  (i)  i 


i=l 


Also, 


a  (u)  =  (l-u)~1{ 


J (w) H' (w) (1-w) dw  +  a. (l-p)H' (p) }, 


and 


2 

a  = 
P 


0 


[a  (u) ]  au. 


It  can  be  shov/n  that 


2  3  , ,  •,  >  ,  2 

°p  =  P  /^-P^p 


-  u 


+  2y  p/A  +  2  J  A^  (-j.)i  +  1/i2i! 


P  P 


i=l 


Then  the  distribution  of  /n((-T/b)  -  y  )/o  approaches  a  standard 

.  P  P 

normal  distribution. 

Also,  k  -»■  -y  as  n  ->  and  the  asymptotic  variance  of 
r ,  n  p 

*  2  2  2 

b  is  b  a  /ny  .  The  asymptotic  variance  of  the  maximum  likelihood 
P  P 

estimator  of  b  is  provided  in  [6]  for  p  =  .  1(.1).9;  and,  since 
this  corresponds  to  the  Cramer-Rao  lower  bound  for  an  unbiased 


estimator  of  b,  the  asymptotic  relative  efficiency  of  b  can  be 

2  *  2 

calculated.  Some  numerical  values  of  k  =  -y  ,  a  ,  n  var  b/b  , 

P  P  P 

and  the  relative  efficiency  of  b  are  presented  in  Table  I. 
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Note  that,  for  complete  samples, 


kn,n  *  xi  -  ln  V!/n 

*  Y  +  E  In  Xn 

=  l  (-1)1  (In  i)  (?) , 
i=l  1 

where  y  denotes  Euler's  constant  and  numerical  values  of  E(ln  X  ) 

'  n 

are  given  in  [18].  Also,  from  [15,  page  71]  k  is  approximated 

11/11 

by  y  +  In  In  n  +  v/ln  n  for  large  n. 


5.  DERIVATION  OF  EXACT  AND  APPROXIMATE  DISTRIBUTIONS 


5.1.  EXACT  DISTRIBUTION  FOR  r  =  2  AND  r  =  3 

The  joint  density  of  (x^,...,xr)  is  given  by  f(x^,...,x  ) 

=  [n!/(n-r)!]{  n  (S/a)  (X^a)  e“1  exp  (-  (X./o)  6]  1  exp  [-  (n-r)  (X  /a)  B] , 
i=l 


0  <  x^  <  ...  <  xr  <  go. 


On  letting  U^  =  X^/Xr,  i  -  l,...,r  -  1,  the  joint  density  of  the 
U^  is  found  to  be 


r-1  R  ,  r-1 

f(u.,...,u  ,)  =  [n!/(n-r)I]  n  8^  xr(r)/[  l  U?  +  n  -  r  +  l]  , 

i=l  i=l 


0  <  u1  <  . . .  <  u  x  <  1, 


*  _  ft  *■  X  A 

which  involves  the  two  quantities  l  U,  and  n  U.  -  exp[-nk  b) 

i=l  1  i*i  1  r,n 

=  S1//B.  Since  is  distributed  independently  of  6,  6  may  be  set 
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equal  to  one  without  loss  of  generality.  For  r  =  2,  S  =  and 
the  distribution  of  S  is  given  directly  (see  also  [7],  [11],  [12]). 

For  r  =  3,  the  distribution  of  S  may  be  determined  by  making  a 
change  of  variable.  The  work  is  simplified  by  noting  first  that 
the  joint  density  of  th  1 s  and  also  the  function  S  are  symmetric  in  the 
variables.  Thus,  for  r  =  3, 


P[S  >  s]  =  P[U^U2  >  s’.-  where  0  <  U.'<  U2  <  1, 


=  iP[U.U_  >  s] ,  where  0  <  U,  <  1 ,  0  <  U„  <  1 , 

Z  X  c  1  4 


Thus, 


1  1 


F3(s)  =  1  -  j 


s/u. 


f  (u^ ,u2) du1du2 , 


Direct  integration  yields  the  result  given  earlier  in  the  paper. 
The  integration  becomes  quite  tedious  for  larger  values  of  r,  so 
that  an  approximation  is  needed. 


5.2.  APPROXIMATE  DISTRIBUTION 

r-i 

The  variable  T  =  -  £  (Y.  -  Y  )/n  takes  on  positive  values, 

i=l  1  r 

and  the  mean  of  nT/b  is  approximately  equal  to  the  variance  of 

nT/b,  especially  for  small  p.  This  is  verified  by  Table  I  for 

2  2 

large  n,  since  var(nT/b)  =  nap,  E(nT/b)  =  nkp,  and  k^  =  ct^  in 
Table  I.  Thus  an  approximate  distribution  with  nearly  the  correct 
first  two  moments  is  obtained  by  assuming  that  2nT/b  is  distribu¬ 
ted  approximately  as  a  chi-square  variable  with  2nk  degrees  of 

r ,  n 

freedom.  The  approximate  percentage  points  were  determined  for 
r  =  2,  n  *  5,10,20;  r  =  3,  n  =  5,10,30,  and  y  =  . 01 , . 05 , . 10 , . 25 , 
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.50, .75, .90, .95  and  .99?  and  the  exact  distributions  were  then 

used  to  determine  the  true  probabilities  for  these  approximate 

percentage  points.  As  seen  in  Table  III,  the  exact  and  approximate 

probabilities  are  in  very  close  agreement. 

This  chi-square  approximation  is  also  consistent  with  the 

2 

asymptotic  results,  at  least  to  the  extent  that  k  =  o  .  This 

P  P 

follows  since  ( (2nT/b)  -  2nk  )/v/4nk_  =  (x  (v)-v)//2v‘,  which 

lr  r*' 

becomes  normally  distributed  as  v  increases?  but 

((2nT/b)  -  2nk  )//4nk  =  /n ( ( T/b )  -  k  )/o  which  corresponds  to 

P  P  P  P 

the  asymptotic  result. 

Thus  the  chi-square  approximation  seems  appropriate  if 
substantial  censoring  is  involved.  Further  work  is  needed  to 
determine  the  amount  of  error  if  r/n  is  near  1. 


TABLE  II 


Comparison  of  b  with  the  BLUE  and  BLIE 


TABLE  III 


Exact  Pr[2nT/b  <  x?  (2nk  )] 

a1~q  r ,n 


JL 

L _ 2 _ 1 

1  3 

n 

5 

10 

20 

5 

10 

30 

01 

.0098 

.0097 

.0099 

05 

.0501 

.0502 

.0500 

10 

.1012 

.1017 

.1003 

25 

.2541 

.2520 

.2510 

.2538 

.2518 

.2508 

50 

.5007 

.  5004 

.  5039 

.5007 

75 

.7417 

.  7454 

.7475 

.  7469 

.7483 

.  7495 

.  8864 

.  8927 

.8962 

.8926 

.  8956 

.  8986 

95 

.9374 

.9433 

.9465 

.9426 

.9460 

.9486 

99 

.9833 

.9865 

.9882 

.9859 

.9878 

.9893 
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CHAPTER  II 


RESULTS  FOR  ONE  OR  MORE  INDEPENDENT  SAMPLES 

1.  INTRODUCTION 

The  results  presented  here  represent  a  continuation  of  the 
previous  chapter  with  particular  application  to  a  problem  con¬ 
sidered  by  Jaech  [7]. 

Jaech  [7]  develops  point  and  interval  estimation  procedures 
for  the  shape  parameter,  g,  if  no  more  than  two  failures  occur 
per  lot-  A  means  for  combining  results  from  two  or  more  lots  is 
also  provided  if  the  shape  parameters  are  assumed  equal.  This 
problem  would  be  of  interest,  for  example,  if  groups  of  items  are 
currently  in  service  for  which  high  reliability  is  required.  Thus, 
as  soon  as  a  few  failures  occur,  the  possible  necessity  of  recal¬ 
ling  all  items  for  replacement  of  degradable  components  would  have 
to  be  considered. 

Procedures  for  point  and  interval  estimation  of  g  based  on 
two  or  more  failures  per  lot  are  presented  in  the  following  sec¬ 
tion  with  a  method  for  combining  results  from  two  or  more  lots. 
Point  estimation  for  a  and  R  is  considered  in  section  3  and  a  con¬ 
servative  lower  limit  for  R  is  given  in  section  4.  Exact  interval 
estimation  procedures  for  a  and  P.  based  on  two  failures  per  lot 
are  derived  in  section  5.  In  section  5  the  relative  expected 
experiment  time  required  to  obtain  a  certain  precision  with  cen¬ 
sored  sampling  as  compared  to  complete  sampling  is  studied.  A 
numerical  example  is  given  in  section  7. 
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2.  ESTIMATION  OF  6  WITH  TWO  OR  MORE  FAILURES 
In  [I]  (Chapter  1) an  estimator  of  b  =  1/8  is  given  as  b  = 

r-1 

-  I  (In  x.  -  In  x  )/nk  =  T/nk  ,  where  the  first  r  failures 

r  t  n  r ,  n 

from  n  items  are  observed  and  the  constants  k  can  be  obtained 

r  ,n 

from  [I]  or  [11] • (The  subscripts  on  the  constants  will  be  suppressed 
hereafter  if  the  meaning  is  clear.)  Also  2T/b  is  distributed  ap¬ 
proximately  as  a  chi-square  variable  with  2nk  degrees  of  freedom, 
at  least  for  r/n  less  than  .5  or  so.  For  r  =  2  the  exact  distri¬ 
bution  of  T/b  is  given  by 

F (t)  =l-net/(et+n-l) 

which  is  approximately  the  exponential  distribution,  1  -  e  t,  for 
large  n.  Thus, 


nk2_n  =  E«/b) 

.00 

=  [1  -  F ( t ) ]  dt 

J0 

=  n  In  ( (n  -  l)/n)  =  1 , 

so  that  1/b  =  1/T ,  which  is  the  estimator  suggested  by  Jaech  [7] 

for  r  =  2.  In  this  case  b  is  the  best  linear  unbiased  estimator 

of  b,  but  as  for  an  exponential  variable,  E(l/b)  is  infinite. 

This  indicates  that  occasionally  very  unreasonable  estimates  will 

occur  if  they  are  based  only  on  a  single  lot  with  r  =  2.  A 

median  unbiased  estimator,  =  c/T,  can  be  found  by  solving 

?[ c/T  <  B]  =  j.  This  gives  c  =  In  [ (2n  -  l)/(n  -  1)]  =  In  2. 

For  r  >  2,Bu  =  (nk  -  l)/(nkb)  is  approximately  an  unbiased  esti- 

2 

mator  of  S,  since  if  U  ~  x  (v)  then  E(l/U)  =  l/(v  -  2)  for  v  >  2. 
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Now  suppose  results  are  available  from  two  or  more  lots  for 
which  a  common  value  of  6  is  assumed.  Suppose  there  are  r ^  ob¬ 
served  failures  from  n ^  elements  in  lot  j,  j  =  l,...,m.  Now  the 

variance  of  a  linear  combination,  £  a.x.  ,  subject  to  [  a.  =  1  is 

j  »  3 

2  2  a  « 
minimized  by  choosing  =  (l/c^  )/£  (l/o^  )  .  Since  Var  (b/b)  == 

2  2 

4nk/4n  k  =  1/nk,  this  leads  to 

„  m  „  m  mm 


b  =  l 


y  n.k.b./)  n.k.  =  V  T ./  \  n.k, 
j=l  3  13  iii  3  3  ^  /j tx  3  3 


m 

as  the  combined  unbiased  estimator  of  b.  Also,  2  £ 

j-1 

tributed  approximately  as  a  chi-square  variable  with 

degrees  of  freedom;  thus  confidence  limits  or  tests 
or  6  are  immediately  available  based  on  one  or  more 


Tj/b  is  dis- 
m 

2  J"  n.k. 
j-1  3  3 

concerning  b 
lots. 


3.  POINT  ESTIMATION  FOR  a  AND  R 


A  closed  form  estimator  of  a  is  given  by  a  = 

{{  l  x.L  +  (n-r  +  l)x  " )  /r  ]  .  This  is  in  the  form  of  the 

i=l  1  r 

usual  maximum  likelihood  estimator  [3]  of  a  with  the  maximum  like¬ 
lihood  estimator  of  3  being  replaced  by  the  simple  closed  form  es¬ 
timator  being  considered  in  this  paper.  Similarly  a  closed  form 

A  ^  O 

estimator  of  reliability  is  given  by  R  =  exp  [  —  ( t/ot )  ]  .  The 
results  of  a  Monte  Carlo  study  of  the  means  and  variances  of 
these  estimators  are  presented  in  Table  1.  The  tabled  values  are 
based  on  2,000  samples  generated  from  a  Weibull  distribution  with 
a  =  1  and  3=1.  Corresponding  values  for  maximum  likelihood  es¬ 
timators  are  available  for  some  cases  [2,5]  and  these  are  included 
for  comparison  purposes.  The  results  of  course  are  applicable 
for  other  values  of  the  parameters  to  the  extent  that  for  both 
methods  of  estimation  3/3  and  (a/a)b  are  pivotal  quantities  with 
distributions  independent  of  both  parameters.  Except  for  the 
simple  estimators  b  and  Bu  there  appears  to  be  substantial  bias 
for  small  r  with  both  methods  of  estimation.  An  unbiased  estima¬ 
tor  of  u  *  In  a  could  be  obtained  for  a  given  n  and  r  by  use  of 
Monte  Carlo  work  if  this  were  deemed  to  be  worthwhile.  For  ex¬ 
ample,  E[B  In  (a/a)]  =  E(ln  a^)  where  denotes  the  estimator 
calculated  from  samples  generated  with  a  =  1  and  6=1.  Thus 
In  a  -  b  E  (In  a^)  is  an  unbiased  estimator  of  In  a. 

To  obtain  a  combined  estimator  from  two  or  more  lots  consider 

~  r — 1 

the  following.  Let  5  =  a6  and  5(B)  =  (  l  x.6  +  (n  -  r  +  l)x  6]/r. 

i=l  1  r 

.  A  2 

It  is  well  known  that  2r£/5  x  (2r)  .  Thus  if  3  were  known, 
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m  „  m 

l  r.5./  l  r.  would  be  the  appropriate  linear  combination  of  the 
j=l  ]  3  j=l  1 

estimators  to  use  to  obtain  a  combined  estimator  of  5  with  minimum 

A 

A  q  A  A  A 

variance.  Since  ct  =  5(6),  this  suggests  using  5  = 

A  C  A 

m  „  6.  m  .  .  b 

T  r.a.  -y  y  r .  as  the  combined  estimator  of  5.  Also  a  =  5 
j=l  D  3  j=l  3  c  c 

* 

and  Rc  =  exp  [- (t/ac)  ] . 
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•?.  INTERVAL  ESTIMATION  OF  R 


Joint  confidence  intervals  for  b  and  R  and  conservative 
lower  limits  for  R  are  given  in  [I]  for  a  single  lot.  Similar 

results  can  be  obtained  for  combined  lots.  As  mentioned  earlier, 

m  *  „  m  m 

2  l  £,($)/£  f\, j  x  i2  l  r.)  and  2  l  T./b  is  distributed  approxi- 
j=l  3  j=l  3  j=l  3 


mately  as 


Also  b .  and 


Cj  (6) 


are  independent  so 


[2£r.?.(B)/C  <  X62(2lr.),Xl?6  (2^.^)  <  2[Tj/b  <  x^U^k.)] 
12  3 


=  (1  -  61)  (1  -  &2  -  63)  , 


2  2 

where  Pfx  (v)  >  X(5  (v)]  =  <5*  This  gives  the  joint  confidence 
region  {R  >  R(b)  ,  b  >  b  >  b} ,  where 

R(b)  =  exp  ( -t 6  x  x; 2  ( 2  l  r.)/2  ]»  r.S.(B)}, 

0  1  J  J  J 

b  =  2  I  T./x62(2  I  n.k.), 
b  =  2  I  T./Xl262  (2  J  n.k.). 

A  conservative  (1  -  6^) (1  -  ~  ^3)  lower  confidence  limit  for  R 

is  given  by  R  =  min_  R(b)  .  As  for  a  single  lot,  R(b)  is  either 
b<b<b 


a  monotonic  function  of  b  or  has  a  single  minimum.  In  particular 

m 

if  t  is  sufficiently  small  and  at  least  if  £  r .  In  (t/x  )  < 
in  j-1  3  j 

-  1  T.,  then  R(b)  is  a  monotonically  decreasing  function  of  b  and 

j=l  3 

R  =  R (E) .  Also  if  t  >  max  xr  ,  then  R  =  R(b) . 

j  j 

It  is  theoretically  possible  to  determine  exact  confidence 
limits  for  R  based  on  R  and  exact  confidence  limits  for  a  based  on 
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A  8 

the  pivotal  quantity  (a/a)  ;  however,  it  has  not  been  possible  in 
general  to  determine  the  necessary  distributions  mathematically. 
Some  Monte  Carlo  work  has  been  carried  out  for  these  quantities 
using  the  less  convenient  maximum  likelihood  estimators  [2],  and 
tables  are  given  in  [8]  for  obtaining  confidence  limits  for  u  or 
b  and  tolerance  bounds  for  the  distribution  bused  on  the  best 
linear  invariant  estimators  of  u  and  b  for  n  =  3(1)25,  r  =  3(l)n. 
Exact  results  for  the  case  r  =  2  are  considered  in  the  next 
section . 
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5.  INTERVAL  ESTIMATION  FOR  a  AND  R 
BASED  ON  TWO  FAILURES  PER  LOT 


For  r  =  2  it  can  be  shown  that  R  is  a  monotonically  increas¬ 
ing  function  of 

v  =  In  (x2/t)/ln  (x,,/x1)  =  In  [  (x2/t)  3]/ [In  (x2/t) L  -  In  (x^/t)^], 

the  distribution  of  which  depends  only  on  6  =  (a/t)^  =  -1/ln  R. 
Letting  Z^  =  (X^/t)0,  then  we  have 

f(21,z2)  =  (n(n  -  l)/02]  exp  [-(z^/d)  -  (n  -  l)z2/e], 

0  <  z^  <  z2  <  ®. 

Making  the  transformation  Y  =  Z^/Z2,  Y2  =  Z2  leads  to 

f(yi,y2)  =  (y2/,02)  exp  t_(YiY2//0)  "  (n  ~  i)y2/eJ» 

0  <  y^  <  1,  0  <  y2  <  ®. 

Thus, 


F  (v )  =  P  [-In  Y2/ln  <  v] 


=  P[Y2  <  Y^]. 


For  vs  0,  F. (0)  = 


_  f1 


rl 


f(y1,y2)  dyx  dy2 


,  -  (n-1) 0  ,  ,  , >  -n/e 

1  -  ne  +  (n  -  l;e 


For  v  <  o,  Fq (v) 


•1/v 


f(y1»y2)  dyi  dy2 


=  F  (0)  -  f  nl-1-.~-il0"(n~1)y2/e  (1  - 


e-*r1/V^,  ay. 


For  v  >  0,  Fq (v) 


re<0>  + 


n(n  -  1) _-(n-I)y„/e ,, 
- § - e  2  (1 


e  y2  'e)  dy. 
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Also,  Fq  (1)  =  ]  -  exp  (-n/e). 

The  above  integrals  were  evaluated  by  numerical  integration 

and  F  (v )  is  tabulated  in  Table  2.  Since  FfJ  (v)  depends  primarily 

u  y 

on  (n  -  1 ) / 0  rather  than  on  n  and  6  separately,  F0 (v)  was  tabu¬ 
lated  as  a  function  of  (n  -  l)/8  and  v  for  n  =  10,  20,  and  100, 
and  linear  harmonic  interpolation  should  be  accurate  for  other 
values  of  n.  For  an  observed  sample  value  v,  a  lower  1-6  con¬ 
fidence  limit,  9_,  for  0  is  the  value  of  0  such  that  F0  (v)  =1-6, 
also  R  =  exp  (-1/Q) .  For  example  suppose  a  90%  lower  confidence 
limit  for  0  is  desired  and  for  n  =  50,  v  =  .3  is  observed.  From 
Table  2  (n  -  l)/0_  =  3  and  0_  =  16.3,  also  R  =  .94.  Of  course  with 
only  two  observed  failures  these  methods  are  not  likely  to  be 
precise  enough  to  be  of  practical  value  unless  results  from  more 
than  one  lot  are  available. 

Although  no  special  techniques  have  been  developed  for  com¬ 
bining  lots  in  this  case,  standard  general  procedures  can  be  used. 

2 

For  example  -2  In  F0 (V)  x  (2)  so  for  m  independent  lots, 
m  _ 

-2  l  In  F. (V. )  xZ (2m) .  For  observed  values  v1f...,v  ,  the 
.  0  l  j.  m 

lower  1-6  confidence  limit  for  9  is  the  value  6_  which  satisfies 

-2  l  In  F0(vi;  =  Xj^Um),  where  P[x2(v)  >  xfi2(v)]  =  6.  The 
i=l  — . 

hypothesis  testing  format  is  somewhat  more  convenient  in  this 

case.  One  would  reject  H  s  9  <  6  in  favor  of  the  alternative 
„  o  -  o 

K  :  9  >  9  if  -2  l  In  F,  (v. )  <  \  (2m). 

•  a  °  i=l  o  1  X 

The  above  results  can  also  be  used  to  determine  confidence 
limits  for  a  by  letting  t  =  a  and  9=1.  That  is,  F^(v)  is  the 
cumulative  distribution  function  for  V  =  In  (X2/a)/ln  (X2/X^)  .  If 


26 


Fl(vv)  =  Y/  then 

V-V 

P[v  <  Vy ]  =  P[u  >  x2(x1/x2)  ]  =  Y / 

vY 

and  x2(x^/x2)  '  is  a  lower  100y%  confidence  limit  for  m 
bilities  for  this  case  are  included  in  Table  2  under  (n 

n  -  1. 


Proba 
-  l)/o 
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6.  COMPARISON  OF  CENSORED  SAMPLING 
AND  COMPLETE  SAMPLING 

Consider  a  test  of  H  s  6  <  6  against  H:  3  >  at  the  6 
level  of  significance.  A  test  for  this  hypothesis  has  been  devel¬ 
oped  in  [9]  for  complete  samples  based  on  the  maximum  likelihood 
estimator,  say  6.  The  null  hypothesis  is  rejected  if  3  >  6q 
where  a ^  satisfies  P [ S/S  >  =  y  and  is  tabulated  in  [9].  The 

power  of  the  test  for  the  alternative  6  is  given  by 

Cl 

P [ 6/ S  >  (B  /3  )&.  .],  If  we  use  the  simple  censored  sample  estimator 
the  null  hypothesis  is  rejected  if  b  <  (1/B0)c. (v),  where  v  =  2nkr  n 
and  c(v)  denotes  the  chi-square  over  degrees  of  freedom  distribu¬ 
tion  [1,4].  The  power  is  given  by  P[c  <  (B_/B  )c  (»)].  Values 
of  v  are  given  in  Table  3  which  would  provide  the  same  power  for 
the  censored  sampling  test  as  would  be  obtained  by  using  the  maxi¬ 
mum  likelihood  method  based  on  a  complete  sample  of  size  N.  Some 
combinations  of  r  and  n  which  would  yield  these  values  of  v  are 
also  given  and  the  relative  expected  experiment  time,  R.E.T.  = 

E  (X  )/E(X..  ,  is  given  in  each  case  for  B  =  1  and  6  =  2.  Note 

that  the  table  indicates  that  the  comparison  does  not  appear  to 
be  very  sensitive  to  the  level  of  significance  or  the  level  of  the 
power  being  considered. 


7.  NUMERICAL  EXAMPLE 

Harter  and  Moore  [6]  give  a  simulated  sample  of  size  40 
from  a  Weibull  population  with  shape  parameter  2,  scale  parameter 
100  and  location  parameter  10;  and,  they  calculated  maximum  like¬ 
lihood  estimates  based  on  the  smallest  10,  20,  30  and  40  observa¬ 
tions  respectively.  After  subtracting  the  value  of  the  location 
parameter  from  each  observation  we  calculated  numerical  values 
which  are  given  in  Table  4  of  the  following  quantities  for 
r  =  2,  10,  20,  30  and  40: 

v  =  2nk  ,  b,  6  =  1/b,  3„  =  <  (v  -  2)/v)6, 

£(.025)  =  x  975(u)/(vb)  =  C  9?5(v)B,  and 

3 (.025)  =  C  q25 (v) 3 /  where  £{6)  denotes  a 

lower  1-6  confidence  limit  for  3.  The 
lower  limit  based  on  maximum  likelihood 
estimates  will  be  denoted  by  3^(6). 

The  hypothesis  Hq:  S  =  =  1  is  rejected  in  favor  of  the  alterna¬ 

tive  6  >  1  at  the  5  significance  level  if  b3Q  <  C^_f(v).  The 
level  of  significance,  6',  at  which  this  hypothesis  could  have  been 
rejected  is  given  in  Table  4.  The  power  of  a  6  level  test  of  Hq  against 
an  alternative  6_  is  denoted  by  P(6,3  )  =  P[C  <  SC.  .(v)]  in 
Table  4.  In  this  example  the  true  reliability  is  .90  at  t  -  32.46. 

A  conservative  .9025  lower  confidence  limit  for  R,  R  = 

2  r“1  " 

min  R(b) ,  where  R(b)  =  exp  {~x  nc(2r)/2[  J  (x./t)^+ 

b ( . 025) <b<E> ( . 025)  i=l  1 

Q 

(n  -  r  +  1) (xr/t)  ] },  was  determined  for  t  =  32.46.  For  r  = 

20,30  and  40  the  condition  In  t  -  In  x  <  -T/r  holds  so  that 


jv 
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R  =  R(b(*025)).  This  is  also  the  limit  for  r  =  10  although  the 
above  condition  does  not  hold  in  this  case. 

The  maximum  likelihood  estimates  were  included  in  the  table 
for  comparison.  Also  confidence  bounds  for  b  and  R  are  available 
from  [10]  for  complete  sampling  and  from  [III]  (Chapter  III)  or  [2] 
for  25%  and  50%  censoring  and  these  are  included.  Although  con¬ 
fidence  bounds  are  included  for  the  complete  sample  case  for  com¬ 
parison  purposes,  it  should  be  recalled  that  the  efficiency  of  b 
and  the  accuracy  of  the  chi-square  approximation  are  not  as  great 
for  large  r/n.  An  estimate  of  a  lower  bound  for  reliability, 

R(b) ,  would  perhaps  also  be  a  useful  statistic  and  it  is  included 
in  the  table. 


TABLE  1 


Monte  Carlo  study  of  means  and  variances  of  estimators; 
o=l,  0=1.  (Tildes  refer  to  maximum  likelihood 
estimators  and  carets  to  the  simple  estimators). 


n 

20  J 

40 

r 

2 

5 

10 

15 

20 

2 

10 

20 

30 

40 

E(b) 

|_  - 

.986 

.995 

.997 

.998 

1.003 

.980 

1.000 

1.000 

1.002 

1.004 

E(b) 

.31 

.81 

.91 

.95 

.96 

E(8) 

1.31 

1.11 

1.06 

1.04 

1.11 

1.046 

1.026 

1.018 

E  (B ) 

1.07 

1.098 

1.060 

1.036 

E<?u> 

1.01 

1.00 

1.01 

1.01 

| 

.999 

.999 

1.000 

1.005 

K(a) 

1.24 

1.007 

.992 

.995 

1.10 

1.000 

.994 

.998 

E(u) 

-.12 

-.060 

-.043 

-.034 

-.06 

-.032 

-.023 

-.017 

E  (u) 

-.34 

-.11 

-.04 

-.02 

ER( . 75) 

.46 

.69 

.738 

.744 

.744 

.42 

.72 

.743 

.746 

.747 

ER( . 75) 

.756 

.752 

.755 

.753 

ER(.90) 

.73 

.89 

.894 

.894 

.892 

.63 

.89 

.896 

.896 

.896 

ER( . 90 ) 

.902 

.904 

.903 

.901 

ER (.95) 

.87 

.94 

.944 

.944 

.943 

.80 

.95 

.947 

.947 

.946 

ER( . 95) 

.950 

.951 

.951 

.950 

V  (b) 

.93 

.218 

.095 

.057 

.042 

.92 

.102 

.044 

.027 

.022 

V(b) 

.23 

.161 

.081 

.051 

.032 

v(bu) 

.89 

.245 

.098 

.057 

.035 

V  ( 6 ) 

.922 

.145 

.074 

.051 

.162 

.054 

.030 

.023 

V(B) 

.036 

.063 

.033 

.019 

V(o) 

1.780 

.140 

.067 

.057 

.580 

.067 

.033 

.028 

v  (u) 

.615 

.136 

.070 

.059 

.290 

.066 

.034 

.029 

V(u) 

.48 

.131 

.07 

.056 

V(R( . 75) ) 

.133 

.030 

.008 

.0068 

.0072 

.148 

.009 

.0035 

.0033 

.0037 

V(R( . 75) ) 

.0062 

.0040 

.0035 

.0031 

V (R( . 90 ) )  | 

,105 

.007 

.003 

.0029 

.0030 

.145 

.  002 

.0016 

.0015 

.0016 

V (R ( . 90 ) ) 

.0023 

.0016 

.0015 

.0012 

V (R ( . 95) ) 

.052 

.002 

.0015 

.0013 

.0013 

.097 

.001 

.0007 

.0007 

.0007 

V (R ( . 95) ) 

.0009 

.0008 

.0007 

.0005 
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TABLE  2 


Values  of  P[ln  ^/tj/ln  ^/X^)  <  v] ,  e  =  (a/t)^  =  -1/ln  R. 

n  =  10 


.  r.-l 

1 

2 

3 

4 

1 

U) 

O 

O 

.006 

.015 

.024 

.032 

O 

O 

CN 

1 

.009 

.023 

.036 

.047 

O 

O 

1 

.017 

.044 

.069 

.090 

-  9.0 

.018 

.049 

.076 

.100 

-  8.0 

.  02  0 

.054 

.085 

.111 

-  7.0 

.023 

.061 

.096 

.125 

-  6.0 

.  027 

.070 

•  111 

.144 

-  5.0 

.031 

.083 

.130 

.168 

-  4.0 

.038 

.101 

.157 

.204 

-  3.5 

.043 

.113 

.176 

.227 

-  3.0 

.049 

.128 

.200 

.257 

-  2.5 

.057 

.149 

.230 

.296 

-  2.0 

.068 

.177 

.272 

.348 

-  1.5 

.084 

.217 

.332 

.421 

-  1.0 

.111 

.281 

.424 

.531 

1 

O 

kT\ 

.161 

.396 

.579 

.704 

-  0.3 

.196 

.469 

.669 

.796 

-  0.1 

.248 

.566 

.773 

.887 

0.0 

.284 

.622 

.823 

.923 

0.1 

.328 

.678 

.864 

.946 

0.3 

.427 

.766 

.913 

.969 

0.6 

.557 

.842 

.946 

.982 

1.0 

.671 

.892 

.964 

.988 

5 

6 

7 

8 

9 

038 

.043 

.047 

.051 

.055 

056 

.063 

.070 

.075 

.080 

107 

.122 

.134 

.144 

.153 

118 

.134 

.147 

.158 

.168 

132 

.149 

.163 

.176 

.187 

149 

.168 

.184 

.198 

.210 

170 

.192 

.210 

.226 

.239 

199 

.225 

.246 

.263 

.279 

240 

.270 

.295 

.316 

.  334 

268 

.301 

.328 

.350 

.370 

302 

.339 

.368 

.393 

.414 

347 

.387 

.420 

.448 

.471 

406 

.452 

.489 

.519 

.544 

489 

.540 

.581 

.614 

.641 

609 

.666 

.709 

.743 

.769 

787 

.842 

.879 

.905 

.923 

872 

.918 

.945 

.963 

.974 

945 

.974 

.987 

.994 

.997 

967 

.987 

.995 

.998 

.999 

979 

.992 

.997 

.999 

1.000 

989 

.996 

.999 

1.000 

994 

.998 

.999 

996 

.999 

1.000 
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TABLE  2  (continued) 


1 

2 

3 

4 

5 

6 

7 

8 

9 

19 

-30.0 

.006 

.015 

.024 

.032 

.038 

.044 

.048 

.052 

.056 

.086 

-20.0 

.008 

.02  2 

.036 

.047 

.057 

.064 

.071 

.077 

.082 

.121 

h-* 

O 

O 

.016 

.043 

.069 

.091 

.109 

.123 

.136 

.147 

.156 

.218 

-  9.0 

.018 

.048 

.076 

.100 

.120 

.136 

.150 

.161 

.172 

.237 

-  8.0 

.020 

.053 

.085 

.111 

.133 

.151 

.166 

.179 

.190 

.262 

.022 

.060 

.096 

.126 

.150 

.170 

.187 

.201 

.214 

.291 

-  6.0 

.026 

.069 

.110 

.144 

.172 

.195 

.214 

.230 

.244 

.329 

-  5.0 

.030 

.082 

.±29 

.169 

.201 

.  228 

.249 

.268 

.284 

.379 

-  4.0 

.037 

.099 

.157 

.204 

.242 

.273 

.299 

.321 

.339 

.446 

-  3.5 

.  042 

.111 

.175 

.228 

.270 

.304 

.332 

.356 

.376 

.489 

-  3.0 

.047 

.126 

.199 

.257 

.304 

.342 

.373 

.399 

.421 

.541 

-  2.5 

.055 

.146 

.229 

.296 

.  349 

.391 

.425 

.453 

.477 

.605 

-  2.0 

.066 

.173 

.270 

.  34  8 

.408 

.455 

.493 

.524 

.550 

.684 

-  ]  .5 

.082 

.213 

.330 

.421 

.490 

.544 

.585 

.619 

.646 

.780 

-  1.0 

.107 

.276 

.420 

.529 

.609 

.668 

.713 

.747 

.774 

.890 

-  0.5 

.155 

.388 

.572 

.700 

.785 

.842 

.880 

.906 

.925 

.983 

-  0.3 

.189 

I 

.459 

.661 

.790 

.869 

.916 

.945 

.963 

.974 

.998 

-  0.1 

.239 

.553 

.763 

.  881 

.941 

.972 

.986 

.993 

.997 

1.000 

0.0 

.274 

.608 

.812 

.916 

.964 

.985 

.994 

.997 

.999 

0.1 

!  .316 

.662 

.852 

.939 

.976 

.991 

.996 

.999 

1.000 

0.3 

!  .412 

.749 

.902 

.963 

.986 

.995 

.998 

.999 

0.6 

.538 

.826 

.937 

.977 

.992 

.997 

.999 

1.000 

1.0 

.651 

.  878 

.957 

.985 

.995 

.998 

.999 
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TABLE  2  (continued) 


n  =  100 


-30.0 

-20.0 

-10.0 

-  9.0 

-  8.0 

-  7.0 

-  6.0 

-  5.0 

-  4.0 


-  3.0 


-  2.0 


-  1.0 

-  0.5 


15  . 
22  . 
43  . 


8  .1 

3  .2 


.189  . 


252  .288 

.385 

.56 

302  .344 

.452 

.63 

335  .380 

.496 

.68 

376  .425 

.548 

.73 

429  .482 

.612 

.79 

497  .555 

.690 

.86 

589  .651 

.785 

.92 

715  .777 

.893 

.97 

880  .926 

.983 

.99 

945  .974 

.998 

1.00 

985  .997 

1.000 

993  .899 

996  .999 

998  1.000 

999 

TABLE  3 


Comparison  between  maximum  likelihood  test  (complete 
sample)  and  simple  censored  sample  test. 


r 

[■' 

; 

f. 


I 

i 


Y 


f 

r 

I 


i 
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TABLE  4 


Numerical  example,  n=40,  6=2,  a=  100. 
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CHAPTER  III 


RESULTS  FOR  CENSORED  SAMPLING  BASED 
ON  THE  MAXIMUM  LIKELIHOOD  ESTIMATORS 

1.  INTRODUCTION  AND  NOTATION 

In  life  testing  experiments  it  is  a  fairly  common  practice 
to  terminate  the  experiment  before  all  items  have  failed.  The 
Weibull  distribution  is  often  used  as  a  model  for  the  observa¬ 
tions  and  when  a  computer  is  available  maximum  likelihood  esti¬ 
mation  of  the  parameters  is  to  be  recommended.  The  tables  pre¬ 
sented  in  this  paper  enable  one  to  set  confidence  limits  on  the 
parameters  and  the  reliability  based  on  the  maximum  likelihood 
estimates  for  selected  censoring  sample  sizes. 

It  is  also  observed  that,  as  in  the  case  with  no  censoring, 
the  maximum  likelihood  estimator  of  the  reliability  is  very  nearly 
unbiased  and  its  variance  is  near  the  Cram6r-Rao  lower  bound. 
Unbiasing  factors  for  he  maximum  likelihood  estimator  of  the 
shape  parameter  aie  given. 

The  form  of  the  Weibull  distribution  function  considered  in 
this  chanter  is 

F(t;b,c)  •  1  -  exp  (-<t/b)c)  for  t  >  0 
where  b  is  the  scale  parameter  and  c  is  the  shape  parameter  The 
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reliability  at  time  t  is  simply  R(t)  =  exp  (-(t/b)c). 

Let  b  and  c  be  the  maximum  likelihood  estimators  of  b  and  c 
and  let  R(t)  be  the  maximum  likelihood  estimator  of  R(t).  Then 
it  is  known  [1,5]  that  the  distribution  of  c/c  and  c  log  (b/b) 
does  not  depend  upon  b  and  c,  although  it  will,  of  course,  depend 
upon  the  sample  size,  n,  and  the  number  of  observations  before 
censoring,  r.  Thus  for  a  given  n  and  r  these  pivotal  functions 
can  be  used  to  test  hypotheses  about  b  and  c  or  set  confidence 
intervals  on  b  and  c.  The  tables  required  when  there  is  no  cen¬ 
soring  are  given  by  Thoman,  Bain  and  Antle  [4],  and  this  paper 
presents  the  tables  when  either  25%  or  50%  of  the  largest  obser¬ 
vations  are  censored.  Moreover,  it  appears  that  linear  interpo¬ 
lation  should  be  adequate  for  censoring  levels  between  those 
given  in  the  tables . 

It  was  shown  [5]  that  the  distribution  of  R(t)  depends  only 
upon  the  values  of  R(t) ,  n  and  r.  Tables  providing  lower  confidence 
limits  for  R(t)  based  upon  m.l.e.'s  from  complete  samples  are 
given  by  Thoman,  Bain  and  Antle  [5],  and  this  paper  presents  the 
tables  needed  when  either  25%  or  50%  of  the  sample  values  are 
censored  from  above. 

The  values  for  each  n  were  obtained  by  simulation  with  8000 
samples  (of  size  n)  used  for  n  =  40,  60,  80,  100  and  120.  The 
8000  samples  were  run  in  two  sets  of  4000  and  the  critical  values 
for  each  set  of  4000  were  compared.  The  critical  values  for  the 
reliability  tables  for  R(t)  >  .9  usually  differed  by  less  than 
.004,  and  so  we  believe  there  is  little  sampling  error  in  these 
tables.  The  critical  values  for  the  other  tables  differed  somewhat 
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more,  those  for  y's  of  .05,  .1,  .9  and  .95  usually  differed  by 
about  .02. 


2.  Inferences  on  the  Parameters 

2.1  Inferences  on  the  shape  parameter 

The  standardized  function  /n  (c/c  -  E(c/c))  was  considered 
in  this  case  because  of  the  convenience  in  the  use  of  the  asymp¬ 
totic  values  and  for  Getter  interpolation  in  the  table.  Table  1 
gives  percentage  points  for  this  quantity  for  selected  cumulative 
probability  levels.  The  asymptotic  percentage  points  were  ob¬ 
tained  from  the  work  of  Harter  and  Moore  [2]  and  are  also  included 
in  the  table.  It  is  seen  from  the  table  that  the  asymptotic 
values  are  approached  quite  slowly.  We  believe  this  is  due  to 
the  lack  of  symmetry  when  the  samples  are  censored  on  one  side, 
and  it  appears  that  the  asymptotic  values  are  not  very  useful  in 
the  censored  case.  Unbiasing  factors  for  c  are  included  in  Table  1. 

Tests  of  hypotheses  concerning  c  or  confidence  intervals  for 
c  based  on  the  function  /n  (c/c  -  E(c/c.)  can  be  easily  developed 
with  the  aid  of  Table  1.  For  example  a  1  -  a  confidence  interval 
for  c  is  given  by 

(c/  (E  (c/c)  +  z1_a/,2//n),c/{E(c/c)  +  z^2//n)  ]  , 
where  the  z^  and  E(c/c)  are  given  in  Table  1. 

2.2  Inferences  on  the  scale  parameter 

Again  as  an  aid  in  interpolation,  percentage  points  for  the 
expression  /n  c  In  (b/b)  are  given  in  Table  2.  Interpolation 
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should  be  fairly  good,  but  as  was  true  for  the  shape  parameter  it 
appears  that  the  approach  to  the  asymptotic  values  is  quite  slow, 
and  for  one  sided  censored  samples  with  n  less  than  120  the  asymp¬ 
totic  values  should  not  be  used.  In  this  case,  for  example,  a 
100(1  -  a)%  confidence  interval  for  b  is  given  by 
[b  exp  (-u^^^/1^"0)  <  k  exp  (-ua^2//nc)  ] . 

3.  Inferences  on  the  Reliability 

In  many  studies  in  which  the  Weibull  distribution  is  used  as 
a  model,  the  primary  interest  is  in  the  reliability  at  some  time 
t,  R (t) .  It  is  fortunate  that  in  spite  of  skewness,  censoring 
and  other  difficulties,  the  m.l.e.  of  R(t)  for  reasonable  values 
of  R(t)  has  negligible  bias  and  its  variance  is  very  close  to 
the  Cram«5r-Rao  lower  bound  for  the  variance  of  an  unbiased  estima¬ 
tor  of  R ( t) .  This  property  was  noted  in  [5]  for  complete  sampling, 
and  it  also  holds  for  censored  sampling.  The  bias  of  R(t)  is 
given  in  Table  3  and  the  variance  in  Table  4.  A  comparison  of 
the  variances  and  the  Cramdr-Rao  lower  bounds  is  given  in  Table  5. 

Table  6  gives  lower  confidence  limits  on  R(t).  These  are 
read  directly  from  the  table  by  entering  the  value  of  R(t)  ob¬ 
served.  A  need  for  these  tables  to  include  high  reliability 
levels  has  been  communicated  to  the  authors,  and  this  accounts 
for  the  number  of  entries  for  reliabilities  near  1  in  the  tables. 

4.  Example 

Harter  and  Moore  [3]  give  a  simulated  sample  of  size  40  from 
a  Weibull  population  with  shape  parameter  2,  scale  parameter  100 
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and  location  parameter  10;  and,  they  calculated  maximum  likelihood 
estimates  based  on  the  smallest  10,  20,  30  and  40  observations, 
respectively.  This  example  with  the  location  parameter  assumed 
known  may  be  used  to  illustrate  the  use  of  the  tables. 

For  r  =  20,  c  =  2.091  and  b  =  83.8.  The  unbiased  estimate 
of  c  is  ( .  911)  (2 . 091)  =  1.90.  Also  for  example,  a  test  of 
Hq:  c  =  1  against  the  alternative  HA:  c  >  1  corresponds  to  a  test 
of  whether  an  exponential  model  is  appropriate,  or  whether  a 
Weibull  model  with  an  increasing  failure  rate  is  needed.  This 
hypothesis  is  rejected  at  the  .05  level  if  /4TT  (c/1  -  1.098)  >  2.95, 
or  if  c  >  1.56.  Thus  the  hypothesis  is  rejected.  A  90%  confi¬ 
dence  interval  for  c  is  given  by 

[2.091/(1.098  +  2.95//40),  2.091/(1.098  -  2.09//40)] 

=  [1.34,  2.72]. 

A  90%  confidence  interval  for  b  is 

[83.8  exp  (-2.16/(2.091) /40) ,  83.8  exp  (3 . 75/ (2 . 091) /40) ] 

=  [71.17,  111.27]. 

In  this  example  the  true  reliability  at  t  =  32.459  is  .90. 

'  7  091 

The  m.l.e.  for  r  =  20  is  R  =  exp  (-(32.459/83.8)  '  )  =  .871. 

From  Table  6  a  lower  90%  confidence  limit  for  R(t)  is  .80. 


TABLE  1 

Percentage  points,  z  ,  such  that  F{/n(c/c  -  E(c/c))  <  z  }  =  y  ,  T  =  /n(c/c) 


TABLE  2 

Percentage  points,  u  ,  such  that  P(/nc  ln(b/b)  <  u  )  =  y 


n 

-V 

.01 

.05 

.1 

.9 

.95 

.99 

1.00 

-2.58 

-1.82 

-1.41 

1.39 

1.80 

2.62 

40 

.75 

-3.29 

-2.25 

-1.69 

1.39 

1.85 

2.61 

.50 

-6.21 

-3.77 

-2.91 

1.63 

2.16 

2.96 

1.00 

-2.48 

-1.78 

-1.38 

1.77 

2.56 

60 

.75 

-3.22 

-2.16 

-1.68 

MR 

1.84 

2.66 

.50 

-5.37 

-3.56 

-2.69 

!  1.67 

2.18 

3.01 

1.00 

-2.51 

-1.76 

-1.37 

1.37 

1.76 

2.49 

80 

.75 

-3.11 

-2.10 

-1.61 

1.43 

1.85 

2.65 

.50 

-5.14 

-3.45 

-2.62 

1.71 

2.16 

3.08 

1.00 

-2.45 

-1.74 

-1.37 

1.35 

1.73 

100 

.75 

-3.12 

-2.09 

-1.60 

1.44 

1.85 

2.61 

.50 

-4.92 

-3.34 

-2.49 

1.78 

2.26 

3.19 

1.00 

-2.44 

-1.73 

I 

-1.35 

1.35 

1.74 

2.48 

120 

.75 

-3.01 

-2.01 

-1.58 

1.45 

1.86 

2.63 

.50 

-4.50 

-3.17 

-2.44 

1.75 

2.27 

3.13 

1.00 

-2.45 

-1.73 

-1.35 

1.35 

1.73 

2.45 

00 

.75 

-2.69 

-1.90 

-1.48 

1.48 

1.90 

2.69 

.50 

-3.69 

-2.61 

-2.03 

2.03 

|  2.61 

3.69 

TABLE  3  BIAS  IN  R(t)  - 


25% 

Censoring 

50% 

Censoring 

R  ( t ) 

40 

60 

80 

100 

120 

40 

60 

80 

100 

120 

.75 

.005 

.003 

.003 

.001 

.002 

.002 

.000 

.001 

.000 

.000 

.80 

.005 

.003 

.003 

.002 

.002 

.004 

.002 

.002 

.001 

.001 

.85 

.004 

.002 

.003 

.001 

.001 

.005 

.002 

.002 

.001 

.002 

.90 

.003 

.001 

.002 

.001 

.001 

.004 

.002 

.002 

.001 

.001 

.925 

.002 

.000 

.001 

.000 

.001 

.003 

.001 

.002 

.001 

.001 

.95 

.001 

.000 

.001 

.000 

.000 

.001 

.001 

.001 

.001 

.001 

.96 

.000 

-.000 

.000 

-.000 

.000 

.001 

.000 

.001 

.000 

.000 

.97 

-.000 

-.000 

.000 

-.000 

-.000 

.000 

.000 

.000 

.000 

.000 

.98 

-.000 

-.000 

-.000 

-.000 

-.000 

-.000 

-.000 

-.000 

-.000 

-.000 

.99 

-.000 

-.000 

-.000 

-.000 

-.000 

-.000 

-.000 

-.000 

-.000 

-.000 

.995 

-.001 

-.000 

-.000 

-.000 

-.000 

-.001 

-.000 

-.000 

-.000 

-.000 

.996 

-.001 

-.000 

-.000 

-.000 

-.000 

-.001 

-.000 

-.000 

-.000 

-.000 

.997 

-.000 

-.000 

-.000 

-.000 

-.000 

-.001 

-.000 

-.000 

-.000 

-.000 

.998 

-.000 

-.000 

-.000 

-.000 

-.Ou.) 

-.001 

-.000 

-.000 

-.000 

-.000 

999 


000  -.000  -.000  -.000  -.000 


000  -.000  -.000  -.000  -.000 


TABLE  4  VARIANCE  OP  R(t) 


25%  Censoring 

50% 

Censoring 

R(t) 

n  40 

60  80  100 

120 

40 

60 

80 

100  120 

.75 

.0035 

.0023  .0018  .0013 

.0012 

.0040 

.0024 

.0019 

.0014  .0013 

.80 

.0030 

.0019  .0015  .0011 

.0010 

.0032 

.0020 

.0016 

.0012  .0011 

.85 

.0023 

.0015  .0012  .0009 

.0008 

.0025 

.0016 

.0013 

.0010  .0008 

.90 

.0015 

.0010  .0008  .0006 

.  0006 

.0016 

.0011 

.0009 

.0007  .0006 

.925 

.0011 

.0008  .0006  ,0005 

.0004 

.0012 

.0008 

.0007 

.0005  .0005 

.95 

.0007 

.0005  .0004  .0003 

.0003 

.0008 

.0006 

.0004 

.0004  .0003 

.96 

.0005 

.0004  .0003  .0003 

.0002 

.0006 

.0004 

.0004 

.0003  .0003 

.97 

.0004 

.0003  .0003  .0002 

.0002 

.0005 

.0003 

.0003 

.0003  .0002 

.98 

.0003 

.0002  .0002  .0002 

.0002 

.0003 

.0002 

.0002 

.0002  .0002 

.99 

.0002 

.0002  .0001  .0001 

.0001 

.0002 

.0002 

.0002 

.0001  .0001 

.995 

.0001 

.0001  .0001  .0001 

.0001 

.0001 

.0001 

.0001 

.0001  .0001 

.996 

.0001 

.0001  .0001  .0001 

.0001 

.0001 

.0001 

.0C01 

.0001  .0001 

.997 

.0001 

.0001  .0001  .0001 

.0001 

.0001 

.0001 

.0001 

.0001  .0001 

.998 

.0001 

.0001  .0001  .0001 

.0001 

.0001 

.0001 

.0001 

.0001  .0001 

.999 

.0001 

.0001  .0001  .0001 

.0001 

.0001 

.0001 

.0001 

.0001  .0001 

i 

i 

! 

! 

i 
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TABLE  5 

COMPARISON 

_ 

OF 

V(R(t)  ) 

WITH 

CRLB 

VARIANCE 

x  104 

(n 

=  40) 

R(t) 

.75 

.80 

.85 

.90 

.925  . 

95  . 

96 

.97 

.98  . 

99 

>.99 

r  =  30 

V (R (t ) ) 

35 

30 

23 

15 

11 

07 

05 

04 

03 

02 

01 

CRLB 

33 

28 

22 

14 

10 

06 

04 

03 

02 

01 

01 

r  =  20 

V  (R  (t)  ) 

40 

32 

25 

16 

12 

08 

06 

05 

03 

02 

01 

CRLB 


35  29  23  16  12  07  05  04  02  01  01 


TABLE  6  A 


90%  LOWER  CONFIDENCE  LIMITS  ON  R(t) 


25% 

CENSORING 

50%  CENSORING 

R(t) 

n  40 

60 

80 

100 

120 

40 

60 

80 

100 

120 

.70 

.623 

.638 

.641 

.650 

.654 

.616 

.639 

.644 

.652 

.655 

.72 

.641 

.657 

.661 

.669 

.673 

.635 

.658 

.663 

.672 

.674 

.74 

.659 

.676 

.681 

.690 

.693 

.653 

.677 

.683 

.691 

.694 

.76 

.678 

.696 

.  702 

.710 

.713 

.674 

.696 

.703 

.711 

.  71  4 

.78 

.698 

.716 

.723 

.731 

.734 

.694 

.716 

.723 

.732 

.734 

.80 

.718 

.737 

.744 

.752 

.755 

.715 

.736 

.744 

.752 

.755 

.82 

.739 

.758 

.766 

.774 

.776 

.737 

.757 

.765 

.773 

.776 

.84 

.761 

.780 

.789 

.796 

.798 

.759 

.779 

.787 

.795 

.797 

.86 

.783 

.802 

.810 

.818 

.821 

.783 

.801 

.810 

.817 

.819 

.88 

.807 

.826 

.833 

.841 

.843 

.807 

.824 

.832 

.839 

.842 

.90 

.832 

.850 

.857 

.864 

.866 

.832 

.847 

.855 

.862 

.864 

.92 

.858 

.875 

.882 

.888 

.890 

.858 

.872 

.879 

.886 

.888 

.94 

.886 

.901 

.907 

.913 

.914 

.886 

.898 

.904 

.910 

.912 

.95 

.901 

.915 

.920 

.925 

.927 

.901 

.911 

.917 

.922 

.924 

.96 

.917 

.929 

.934 

.939 

.940 

.917 

.925 

.930 

.935 

.937 

.97 

.938 

.943 

.947 

.952 

.953 

.933 

.940 

.944 

.949 

.951 

.98 

.951 

.959 

.963 

.966 

.967 

.951 

.956 

.959 

.964 

.965 

.99 

.971 

.976 

.979 

.981 

.982 

.971 

.974 

.977 

.979 

.980 

.9925 

.977 

.981 

.984 

.986 

.986 

.977 

.979 

.982 

.984 

.985 

.995 

.983 

.987 

.989 

.990 

.990 

.983 

.985 

.987 

.988 

.989 

.996 

.986 

.989 

.990 

.992 

.992 

.986 

.987 

.939 

.990 

.991 

.997 

.989 

.992 

.993 

.994 

.994 

.989 

.989 

.991 

.992 

.993 

.998 

.992 

.534 

.995 

.995 

.996 

.992 

.992 

.994 

.995 

.995 

.9985 

.993 

.995 

.996 

.996 

.997 

.993 

.994 

.995 

.996 

996 

.999 

.994 

.996 

.997 

.998 

.998 

.994 

.995 

.996 

.997 

.997 
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TABLE  6  B 


95%  LOWER  CONFIDENCE  LIMITS  ON  R(t) 


251 

CENSORING 

50% 

CENSORING 

R(t) 

n  40 

60 

80 

100 

120 

40 

60 

80 

100 

120 

.70 

.594 

.626 

.624 

.625 

.643 

.600 

.623 

.  623 

.639 

.646 

.72 

.613 

.644 

.643 

.647 

.662 

.614 

.  641 

.  647 

.659 

.664 

.74 

.632 

.662 

.664 

.669 

.681 

.632 

.660 

.667 

.678 

.683 

.76 

.651 

.680 

.684 

.691 

.701 

.651 

.679 

.686 

.698 

.702 

.78 

.671 

.699 

.705 

.713 

.722 

.671 

.699 

.707 

.719 

.722 

.80 

.692 

.719 

.726 

.736 

.743 

.691 

.719 

.727 

.741 

.742 

.82 

.714 

.740 

.748 

.759 

.764 

.712 

.740 

.  749 

.761 

.762 

.84 

.737 

.761 

.771 

.782 

.786 

.  734 

.761 

.  771 

.782 

.784 

VD 

00 

.760 

.784 

.795 

.  806 

.809 

.757 

.  784 

.  793 

.805 

.806 

.88 

.785 

.808 

.819 

.830 

.  832 

.781 

.  807 

.  817 

.827 

.829 

.90 

.811 

.833 

.844 

.854 

.856 

.807 

.831 

.841 

.851 

.852 

.92 

.839 

.860 

.870 

.879 

.881 

.834 

.857 

.  866 

.  876 

.877 

.94 

.869 

.888 

.897 

.904 

.906 

.  863 

.883 

.  892 

.902 

.903 

.95 

.  885 

.903 

.911 

.917 

.920 

.878 

.879 

.906 

.915 

.917 

.96 

.902 

.019 

.926 

.932 

.  933 

.894 

.913 

.920 

.929 

.931 

.97 

.920 

.935 

.941 

.946 

.948 

.913 

.929 

.936 

.946 

.947 

.98 

.940 

.953 

.957 

.962 

.963 

.933 

.947 

.952 

.960 

.961 

.99 

.964 

.972 

.976 

.978 

.979 

.957 

.968 

.971 

.975 

.977 

.9925 

.370 

.978 

.981 

.  9  §  3 

.984 

.965 

.974 

.977 

.980 

.981 

.995 

.978 

.984 

.986 

.988 

.988 

.973 

.980 

.983 

.  985 

.986 

.996 

.981 

.986 

.988 

.990 

.990 

.976 

.983 

.985 

.988 

.989 

.997 

.985 

.989 

.991 

.992 

.992 

.980 

.986 

.988 

.990 

.991 

.998 

.988 

.992 

.993 

.994 

.995 

.985 

.990 

.992 

.993 

.  994 

.9985 

.991 

.994 

.995 

.996 

.996 

.987 

.992 

.993 

.994 

.995 

.999 

.993 

.995 

.996 

.997 

.997 

.990 

.994 

.995 

.996 

.996 
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TABLE  6C 


98%  LOWER  CONFIDENCE  LIMITS  ON  R(t) 


25% 

CENSORING 

r 

50% 

CENSORING 

R(t) 

n  40 

60 

80 

100 

120 

40 

60 

80 

100 

120 

.  70 

.571 

c  - 

.610 

.622 

.  625 

.590 

.606 

.613 

.628 

.631 

,72 

.589 

.615 

.629 

.641 

.6  45 

.606 

.623 

.632 

.647 

.649 

.  74 

.601 

.635 

.648 

.661 

.666 

.622 

.641 

.650 

.665 

.669 

.76 

.627 

.654 

.  668 

.681 

.687 

.639 

.660 

.67^ 

.685 

.689 

.78 

.647 

.675 

.688 

.702 

.708 

.657 

.679 

.689 

.704 

.708 

.  80 

.669 

.696 

.709 

.723 

.730 

.676 

.698 

.710 

.725 

.729 

.  82 

.690 

.718 

.731 

.745 

.752 

.695 

.719 

.731 

.746 

.754 

.84 

.713 

.741 

.753 

.768 

.775 

.716 

.740 

.752 

.768 

.773 

r\ ,  • 

•  OU 

.737 

.764 

.777 

.791 

.798 

.737 

.763 

.775 

.790 

.792 

.88 

.763 

.789 

.801 

.815 

.  822 

.760 

.  78C 

.799 

.813 

.819 

.90 

.790 

.816 

.827 

.840 

.846 

.735 

.811 

.824 

.837 

.848 

.92 

.819 

.844 

.854 

.866 

.871 

.812 

.838 

.850 

.863 

.871 

.94 

.851 

.873 

.883 

.893 

.896 

.842 

.866 

.877 

.889 

.893 

.95 

.868 

.889 

.898 

.908 

.911 

.858 

.881 

.892 

.903 

.906 

.96 

.886 

.906 

.914 

.923 

.926 

.875 

.897 

.907 

.918 

.922 

.97 

.906 

.924 

.931 

.938 

.941 

.895 

.915 

.924 

.934 

.937 

.98 

.928 

.943 

.950 

.955 

.957 

.917 

.935 

.943 

.951 

.953 

.99 

.954 

.966 

.970 

.974 

.976 

.945 

.959 

.964 

.971 

.972 

.9925 

.962 

.972 

.976 

.980 

.981 

.953 

.966 

.970 

.976 

.978 

.995 

.971 

.979 

.983 

.985 

.986 

.963 

.974 

.977 

.982 

.984 

.996 

.975 

.983 

.985 

.988 

.989 

.967 

.978 

.981 

.985 

.987 

.997 

.979 

.986 

.988 

.990 

.991 

.972 

.982 

.984 

.988 

.989 

.998 

.984 

.990 

.991 

.993 

.994 

.978 

.986 

.988 

.991 

.992 

.9985 

.987 

.992 

.993 

.994 

.995 

.982 

.989 

.990 

.993 

.994 

.999 

.990 

.994 

.995 

.996 

.996 

.986 

.992 

.993 

.995 

.995 
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17.BLE  6D 


95%  LOWER  CONFIDENCE  LIMITS  ON  R(t) 


25%  CENSORING 

_ 

50%  CENSORING 

R(t) 

n  40 

60 

80 

100 

120 

40 

60 

80 

100 

120 

.70 

.555 

.585 

.601 

.613 

.623 

.566 

.  590 

.609 

.613 

.615 

.72 

.574 

.603 

.620 

.636 

.641 

.582 

.607 

.626 

.633 

.636 

.74 

.592 

.622 

.638 

.655 

.661 

.599 

.624 

.643 

.652 

.  656 

.76 

.612 

.642 

.658 

.674 

.680 

.617 

.643 

.661 

.672 

.677 

.78 

.632 

.662 

.678 

.694 

.701 

.636 

.662 

.679 

.693 

.698 

.80 

.652 

.684 

.699 

.715 

.  721 

.655 

.6  81 

.698 

.714 

.720 

.82 

.674 

.705 

.720 

.736 

.743 

.675 

.702 

.718 

.735 

.742 

.84 

.697 

.728 

.743 

.759 

.765 

.696 

.723 

.739 

.758 

.765 

.86 

.722 

.752 

.766 

.782 

.788 

.718 

.746 

.761 

.  780 

.788 

.88 

.747 

.777 

.791 

.806 

.  812 

.742 

.770 

.  784 

.  804 

.812 

.90 

.775 

.804 

.816 

.831 

.837 

.768 

.796 

.  809 

.829 

.836 

.92 

.803 

.832 

.844 

.858 

.863 

.795 

.823 

.836 

.  854 

.861 

.94 

.838 

.863 

.873 

.886 

.890 

.826 

.853 

.865 

.  881 

.886 

.95 

.855 

.879 

.889 

.901 

.905 

.843 

.869 

.881 

.  895 

.899 

.96 

.874 

.896 

.906 

.916 

.920 

.861 

.886 

.897 

.910 

.914 

.97 

.895 

.915 

.923 

.933 

.936 

.881 

.905 

.915 

.927 

.931 

00 

a\ 

.918 

.936 

.943 

.951 

.954 

.904 

.926 

.935 

.945 

.949 

.99 

.947 

.960 

.966 

.971 

.973 

.  934 

.953 

.959 

.965 

.967 

.9925 

.956 

.968 

.973 

.977 

.978 

.943 

.960 

.966 

.971 

.974 

.995 

.966 

.975 

.980 

.984 

.984 

.954 

.969 

.977 

.978 

.982 

.996 

.970 

.979 

.983 

.986 

.987 

.959 

.973 

.  978 

.981 

.985 

.997 

.975 

.983 

.986 

.989 

.989 

.966 

.*77 

.982 

.985 

.988 

.998 

.981 

.987 

.990 

.992 

.992 

.974 

.982 

.986 

.988 

.991 

.9985 

.984 

.990 

.992 

.993 

.994 

.978 

.975 

.988 

.991 

.99  3 

.999 

.987 

.992 

.994 

.995 

.996 

.983 

.988 

.991 

.99  3 

.994 
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